捷联式航空重力测量系统的数据处理方法比较与分析

重力数据处理对获取高精度重力异常值有着重要作用，是重力测量的核心技术。重力仪在搭载运动载体进行重力测量时，载体的高频振动对重力测量数据和GPS数据均会产生不可避免的干扰，导致提取的重力异常粗值含有大量高频噪声。围绕重力数据的处理方法这一核心技术，介绍了FIR低通滤波、零相移滤波、标准Kalman滤波、正反Kalman滤波4种滤波方法的基本原理，运用这4种方法处理了SAG捷联式重力仪的某次实际飞行测量数据，比对了基于SINS/GPS组合导航和SINS/DGPS组合导航的重力测量数据处理结果。通过对本次试验重复测线内符合精度进行比对，验证了4种方法的可行性和优劣性，同时验证了SAG重力仪的测量精度。     

Substitution of satellite-based land surface temperature defective data using GSP method

Land surface temperature (LST) as an important environmental variable provides valuable information for earth environmental system modelling. Currently, LST is obtained through satellite thermal sensors at various spatial and temporal resolutions. Although spatially continuous satellite-based LST measurements are intended to overcome the shortcomings of sparse ground-based LST measurements, LST images often contain anomalous values due to the existence of clouds or sensor malfunctioning. The problem becomes more serious where the users deal with high spatial resolution characterized by low temporal resolution. This study examines the capability of a newly developed graph signal processing (GSP) method using two-dimensional single-date thermal data. For this purpose, four Landsat/TIRS datasets are analyzed. The data of five elliptical regions on thermal images are eliminated and then reconstructed through the GSP method and using the LST values of the enclosing rectangles containing the ellipsoids. The results indicate that the temperature variation determined by the GSP method generally conforms to the original image LST values. According to a correlation test conducted on the original image LST and those obtained through the GSP method, the values vary from 58% to 95%, which is an above-the-average rate (RMSE from 0.69 to 2.27). The statistical analysis of the original image LST in both the elliptical regions and the enclosing rectangles containing the ellipsoids indicates that an increase in the variance of LST data causes an increased error in the calculation of temperature by the GSP method, and vice versa. The results of the analysis of variance (ANOVA) and Duncan test indicated that an increase in the number of the non-zero spectral bins would result in increased RMSE values for all the dates and the regions. Moreover, the model errors were significant at the 0.05 level across all the image date and five elliptical study regions. Based on the results, the use of this method is recommended for the reconstruction of LST missing values, where dissimilarity of atmospheric conditions limits the use of other methods that depend on the time series data of various dates and a great deal of data calculation.     

应用后缘装置的跨声速层流翼型多岛优化设计

高雷诺数状态下,自然层流技术(Natural laminar flow,NLF)是减小机翼表面湍流摩擦阻力的有效方法。然而由于层流翼面上大范围顺压梯度的存在使得后缘处的恢复压差更大,产生更强的激波。因此在减小摩擦阻力的同时又增加了激波阻力。本文采用后缘装置(Trailing edge device,TED)来控制翼型后缘处的激波强度,基于线性稳定性理论(Linear stability theory,LST)的eN方法对流动进行转捩判断,进而应用多岛并行多目标进化算法(multi-objective evolutionary algorithm,MOEA)以获得大范围层流区域和弱化激波强度为目标对翼型进行优化设计。优化结果表明合作均衡策略耦合进化算法可以快速地捕捉到该多目标问题的Pareto阵面解,阵面上翼型的波阻力和摩擦阻力性能较初始翼型大大改善。同时,采用后缘装置控制激波强度时,无论在设计点还是偏离设计点时,优化后翼型均具有良好的升阻力特性和鲁棒性。     

Direction finding for two-dimensional incoherently distributed sources with Hadamard shift invariance in non-uniform orthogonal arrays

This paper proposes a novel algorithm for Two-Dimensional (2D) central Direction-of-Arrival (DOA) estimation of incoherently distributed sources. In particular, an orthogonal array structure consisting of two Non-uniform Linear Arrays (NLAs) is considered. Based on first-order Taylor series approximation, the Generalized Array Manifold (GAM) model can first be established to separate the central DOAs from the original array manifold. Then, the Hadamard rotational invariance relationships inside the GAMs of two NLAs are identified. With the aid of such relationships, the central elevation and azimuth DOAs can be estimated through a search-free polynomial rooting method. Additionally, a simple parameter pairing of the estimated 2D angular parameters is also accomplished via the Hadamard rotational invariance relationship inside the GAM of the whole array. A secondary but important result is a derivation of closed-form expressions of the Cramer-Rao lower bound. The simulation results show that the proposed algorithm can achieve a remarkably higher precision at less complexity increment compared with the existing low-complexity methods, which benefits from the larger array aperture of the NLAs. Moreover, it requires no priori information about the angular distributed function.     

Latitudinal distribution of solar coronal active regions

We studied the cyclic evolution of the latitudinal distribution of solar coronal active regions based on daily images from SOHO EIT for the period 1995–2017. Fully automated software was used, which included the following steps: initial preparation of images in the data series, normalization of histograms and correction of limb brightening, segmentation of images using threshold intensity values obtained from their histograms, scanning of segmented images in heliographic coordinates and obtaining profiles of latitudinal distribution of coronal active regions for each image of the data series. From the output data, we obtained a temporary change in the latitudinal distribution profiles and the migration of activity centers on the solar disk. From the period of minimum activity to the next minimum in both hemispheres, activity centers begin to migrate from high latitudes towards the equator. At the same time, the general center of activity repeatedly changes the direction of migration. The latitudinal distribution of the so-called presence factor of coronal active regions closely resembles the magnetic butterfly diagram, which proves their direct causal relationships. Variations in the presence factor of coronal active regions are correlated with cyclic variations in the sunspot daily numbers.     

Real-time clock comparison and monitoring with multi-GNSS precise point positioning: GPS,GLONASS and Galileo

The precise point positioning (PPP) technique is widely used in time and frequency applications. Because of the real-time service (RTS) project of the International GNSS Service, we can use the PPP technique for real-time clock comparison and monitoring. As a participant in the RTS, the Centre National d’Etudes Spatiales (CNES) implements the PPPWIZARD (Precise Point Positioning with Integer and Zero-difference Ambiguity Resolution Demonstrator) project to validate carrier phase ambiguity resolution. Unlike the Integer-PPP (IPPP) of the CNES, fixing ambiguities in the post-processing mode, the PPPWIZARD operates in the real-time mode, which is also called real-time IPPP (RT-IPPP). This paper focuses on applying the RT-IPPP for real-time clock comparison and monitoring. We review the principle of real-time clock comparison and monitoring, and introduce the methodology of the RT-IPPP technique. The observations of GPS, GLONASS and Galileo were processed for the experiments. Five processing modes were provided in the experiment to analyze the benefits of ambiguity resolution and multi-GNSS. In the clock comparison experiment, the average reduction ratios of standard deviations with respect to the G PPP mode range from 9.7% to 35.0%. In the clock monitoring experiment, G PPP mode can detect clock jumps whose magnitudes are larger than 0.9 ns. The RT-IPPP technique with GRE PPP AR (G) mode allows for the detection of any clock jumps larger than 0.6 ns. For frequency monitoring, G PPP mode allows detection of frequency changes larger than 1.1 × 10⁻¹⁴. When the RT-IPPP technique is applied, monitoring with GRE PPP AR (G) mode can detect frequency changes larger than 6.1 × 10⁻¹⁵.     

平台式重力仪航空重力测量数据后处理技术研究

针对航空重力测量对数据后处理的要求,结合平台式重力仪系统的组成与工作原理,在重力异常信号的频域特征基础上,设计出一套有严格频域规范的航空重力数据后处理流程,并对某海域实测数据进行了处理。结果表明,重力异常值内符合精度为12条东西测线达到0.43mGal/140s和0.84mGal/100s,4条南北测线达到0.39mGal/140s和0.79mGal/100s;以同测线GT-2A单次测量结果为标准值统计标准差,12条东西测线为0.72mGal/140s和0.98mGal/100s,4条南北测线为1.41mGal/140s和1.53mGal/100s。结果表明,设计的航空重力测量数据后处理方法可以实现重力异常信号的高精度提取。     

基于STM32的视频流图像更新方法的研究

对基于STM32系统低速(工作主频72MHz)视频流处理显示更新方案的研究。结合图像处理基本理论,分析传统的显示方法的不足,并提出一种新的视频流图像更新的方法,通过一定的实验数据和理论结果,验证该方法可行性。     

改进层次聚类法在复合材料敲击检测中的应用

伴随复合材料的迅速发展,其在航空航天方面的应用也日益普遍,但因复合材料的各向异性,在制造过程中会导致某些缺陷的产生。敲击检测作为一种实时原位的无损检测,在复合材料的检测中效果明显。由于敲击检测尚未找到合理的数据处理方法,所以其应用范围受到限制。从解决敲击检测数据处理方法的角度出发,提出将改进层次聚类法应用于敲击检测中,并在应用实例环节对该方法进行了应用,通过实验结果对比,表明改进层次聚类方法可以解决敲击检测数据处理的问题。